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Abstract. We present a FUSE spectroscopic study of the nearby gas-rich, metal-deficient blue compact dwarf 
(BCD) galaxy IZw 36. Atomic hydrogen and many metal lines are observed in absorption against the stellar 
continuum of young, massive stars embedded in the ionized region. Profile fitting of absorption lines allowed us 
to determine abundances and investigate the chemical composition of the neutral gas. This study present strong 
evidences that nitrogen is 16^g times less abundant in the neutral gas with respect to the ionized gas (all uncer- 
tainties are 2a). Similarly, the oxygen abundance estimated using phosphorus as a tracer is lower in the neutral 
gas by a factor 8I5 7 . We also find that argon is underabundant by a factor 32l? and that log (Ar i/O l)<— 3.0 
which is inconsistent with the solar value — 2.1±0.1 (Lodders 2003), implying that argon is likely ionized into Aril 
in the neutral medium. 
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1. Introduction 

Blue compact dwarfs (BCDs) are gas-rich and relatively 
■ unevolved objects, as suggested by their low metallic- 
ity (between ~ Zq/50 and ~ Zq/3). Their star forma- 
tion history is thought to be dominated by short, in- 
tense bursts separated by long quiescent periods (Searle 
& Sargent 1972). The present starburst episode gives 
an opportunity to compare the chemical abundances of 
the H 11 region where stars recently formed with those 
of their surrounding extended neutral gas. One of the 
issues is to know whether metal abundances as mea- 
sured from H 11 regions reflect the true metallicity of star- 
forming galaxies and whether the neutral phase is less 
chemically evolved than the ionized phase. Such an at- 
tempt has been done in the past using Hubble Space 
Telescope (HST), with inconclusive results because of 
the limited wavelength coverage and the difficulties in 
analyzing the saturated A1302.2 A Oi line (Kunth et 
al. 1994 and Pettini & Lipman 1995). The Far Ultraviolet 
Spectroscopic Explorer (FUSE, Moos et al. 2000) gives 
access to many transitions of critical species such as H 1, 
Ni, Oi, Sill, Pll, Arl, or Fell, allowing one to revisit the 
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chemical composition in neutral regions and the metal en- 
richment processes in nearby galaxies. 

The fate of the metals newly produced by starbursts 
in BCDs is not settled yet. Kunth & Sargent (1986) have 
emitted the hypothesis that the H11 region of IZw 18 has 
enriched itself over the timescale of the burst (i.e. few 
10 6 yr) . The H 11 region could be self-polluted, resulting in 
larger abundances than the surrounding neutral medium. 
Another possibility is that, once released by massive stars, 
metals remain in a hot phase, being unobservable immedi- 
ately through optical and UV emission lines. They would 
then cool into molecular droplets which later on settle 
onto the disk (Tenorio-Tagle 1996). Mixing with the in- 
terstellar medium (ISM) has to await photodissociation 
of the droplets by the UV radiation of another genera- 
tion of stars before their material changes the interstellar 
medium composition. 

The study of various elements in both ionized and 
neutral gas in star-forming galaxies should help to con- 
strain the enrichment processes and identify parameters 
that lead to possible abundance differences such as deple- 
tion on dust grains, ionization state correction and unpro- 
cessed (less chemically evolved) gas in the line of sight. 

In this paper we investigate the ISM of IZw 36 through 
absorption lines of species standing in the line of sight of 
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the ionizing stellar complex within the BCD. IZw 36 is 
a metal-deficient BCD with Z=Z©/14 in the Hn region 
(Viallefond & Thuan 1983, hereafter V83). The Hi distri- 
bution shows a core-halo structure. The ionizing cluster 
is associated with the H I core which contains half of the 
total H I mass (V83) . The H I halo is diffuse and contains 
several clumps. 

In order to compare the chemical composition of the 
neutral gas with respect to the ionized gas, we will con- 
sider the abundances of nitrogen, oxygen, silicon, argon, 
phosphorus, and iron. The observations are described in 
Sect. and the data analysis is explained in Sect. OH The 
composition of the neutral gas of IZw 36 is derived in 
Sect. ~Q In Sect. we discuss the abundance differences 
between the ionized and the neutral gas. 

2. Observations 

IZw 36 was observed with FUSE on 2001, January 13 (ob- 
servation Q2240101) with an integration time of 19350 sec- 
onds (7 sub-exposures) and on 2000, May 9 (observation 
PI 072201) with 6290 seconds (3 sub-exposures). Because 
of the angular extent of the ionizing cluster, the largest 
aperture LWRS was chosen (30" x 30" ) . At a distance of 
6 Mpc (Tabled, 1" corresponds to 30 pc. 

Table 1. Characteristics of IZw 36. 



Parameter 


Value 


a(2000), 5(2000) 


12 h 26 m 16 s , +48°29'37" 


Mb 


-14.07 


v a (km s^ 1 ) 


281±4 


D b (Mpc) 


5.8 < D < 7.9 


Z c 


« Z /14 



a Heliocentric radial velocity. 

b Distance (Schulte-Ladbeck et al. 2001). 

c Metallicity (V83). 



Data were recorded through the two LiF channels 
(^1000-1200 A) and the two SiC channels («900-1100 A) 
and processed by the pipeline Calf use 2.2.0. The spec- 
trum of each channel results from the co-addition of sev- 
eral individual sub-exposures. The typical signal-to-noise 
ratio (S/N) per resolution element of the co-added spectra 
ranges from 2 below 1000 A up to 4 above 1000 A. The 
final spectrum shows two absorption lines systems sepa- 
rated by approximately 1 A. The first at ss— 50 km s _1 
is due to the local ISM of the Milky Way. The second 
is at 270^| km s _1 , which matches with the heliocentric 
velocity of IZw 36 (Tabled). 

3. Data analysis 

To analyze the data, we used the profile fitting procedure 
Owens developed at the Institut d'Astrophysique de Paris 



by Martin Lemoinc and the FUSE French team. This pro- 
gram returns the most likely values of many free param- 
eters such as temperatures, Doppler widths, velocities, or 
column densities through a minimization of the difference 
between the observed and computed profiles of absorption 
lines. Furthermore, the version we used allows changes of 
the background level, the continuum, and the line broad- 
ening. The continua were fitted by zero- to fourth-order 
polynomials, depending on the spectral region. All the co- 
efficients of the polynomials were free. An example of fit- 
ted profiles can be seen in Fig. Q] 

We consider independently several groups of species, 
each group being defined by its turbulent velocity, its 
temperature, and its heliocentric velocity. As shown in 
Table El one group defines species supposed to be mainly 
present in the neutral gas (see also Sect. I4.2JI . The two 
other groups contain respectively species of the ionized 
phase and the hot gas (giving the absorption lines of O Vl). 
Molecular hydrogen is discussed separately in Sect. 14. ll 

3.1. The line broadening 

Taking the observed total line broadening otot as a free 
parameter, we used Owens to find its most likely value 
in IZw 36 spectra, supposing that a tot is constant in 
the FUSE spectral range. We find er fot =27±3 pixels 
(»50 km s" 1 at A = 1000 A). 

The total broadening of the unsaturated absorption 
lines has several origins: 

— The instrumental line spread function. We could not 
use the H2 lines from IZw 36 and the Galactic compo- 
nent, which usually allow a good estimation of the in- 
strumental broadening, since they are not detected in our 
spectra (Sect, l4~l.ll) . For a bright point -like source, the 
full width at half maximum (FWHM) is tJi ns t — H pixels 
(Hebrard et al. 2002). 

— The misalignments of individual sub-exposures. The fi- 
nal spectrum of each detection channel is obtained by co- 
additions requiring wavelength shifts of individual sub- 
exposures. This unavoidably introduces some misalign- 
ments because of the low S/N ratio of each exposure. 
We estimate that such misalignments cause an addition- 
nal broadening of the LSF of a a dd ~10 pixels. 

— The spatial distribution of the UV bright stars within 
the slit. HST/FOC images (Deharveng et al. 1994) reveal 
that the main concentration has an extent of «4" , corre- 
sponding to a wavelength smearing of o~ spa t ~7 pixels. 

— The main source of line broadening is the velocity dis- 
tribution of the absorbing clouds standing in the multiple 
lines of sight, whose width is a priori unknown but can be 
deduced for an unsaturated line by <T c i ou ds = {otot ~ a ?nst ~ 
a 2 spat ) 1 / 2 =22 pixels i.e. 40 km s _1 . 

The influence of the resulting total broadening a tot 
on the determination of column densities is discussed in 
Sect-UZO 
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Table 2. List of the lines detected in IZw 36. MW is for the Milky Way, A is the rest wavelength. Velocities of each species 
have been computed independently; typical errors are ±10km s _1 . Species are distributed into three groups (see text), b and v 
are calculated using line profiles of all species within a group. 



Species 


Velocity 


A 


comment 




(km s" 1 ) 


(A) 








First group = neutral region 


of IZw 36. 6=12.2+^ km s~~\ u=270±g km s~ x 


Hi 




1025.72 


Ly/3, presence of damping wings (Fig. ^i) 






937.80, 930.75, 926.23, 923.15 


resp. Lye, Ly£, Lyr/, Ly8 


Cii 


259 


1036.34 


strongly saturated, blended with C n* from MW and VI from MW 


Ni 


265 


1134.98 


not saturated (Fig.|SJ 






963.99 


not saturated, blended with N I from MW 








and P II from IZw 36 






954.10, 953.97, 953.65, 953.42 


not saturated, blended with each other 


Oi 


272 


1039.23 


strongly saturated 






976.45 


saturated, blended with C in from MW 






950.89 


saturated, blended with P IV from IZw 36 


Sin 


261 


1020.70 


slightly saturated 


Pn 


270 


1152.82 


detected 






963.83 


blended with N I from MW and N I from IZw 36 


An 




1048.22 


not detected (Fig.|3Jl 


Fen 


270 


1144.94 


not saturated (Fig.^j) 






1143.23, 1142.37, 1125.45 


not saturated 






1121.97, 1096.88, 1063.18 


not saturated 




Second 


group: species with ionization potentials>29 eV. 6=24. 3lig % km s 1 , t;=260^* 2 km s 1 


Cm 


261 


977.02 


detected 


Nil 


254 


1083.99 


barely detected 


Piv 


254 


950.66 


blended with O I from IZw 36 


Sin 


248 


1012.50 


detected 


Fein 


264 


1122.53 


detected 


Third group = hot gas. 6=51.4+^;^ km s"\ «=235±^ km s" 1 


Ovi 


235 


1037.62, 1031.93 


detected 



4. Composition of the neutral gas 

4.1. Column densities 

Column densities are given in Table [21 Althrough the 
paper, error bars are 2a. The multiple lines of sight to 
the massive stars that contribute to the UV continuum 
may cross clouds with different chemical composition, so 
that the reported column densities represent global values 
among all the clouds. 

Errors on the column densities essentially depend on 
where the lines are located on the curve of growth. In par- 
ticular, saturated lines can give column densities spanning 
3 orders of magnitude when b is not constrained by lines 
of other species standing in the same gas phase. Errors 
are computed using the A^ 2 method: x 2 is calculated 
as a function of the column density. We use the differ- 
ence A% 2 =x 2 — Xmim where Xmin i s the best fit to derive 
the error bars (2a error bars correspond to A% 2 =4). This 
method is also used to evaluate error bars for b and the line 
broadening. Figure [21 shows an example of a A\ 2 curve. 

The determination of the H I column density is not 
affected by the line broadening a tot- However, we have 



3.2. The turbulent velocity 



The turbulent velocity parameter b of the neu- 
tral gas is well constrained by both saturated- 
profile and Doppler-profile lines. From profile fit- 
ting we obtain b=12.2^\ km s _1 (2a error bars), 
corresponding to a velocity dispersion FWHM of 
Orient™; =2 Vhi 2 6=20.3 km s _1 . This result is consis- 
tent with the velocity dispersion measured from 21 cm 
radio observations (021 cm «45 km s _1 , V83) and the de- 
rived velocity distribution of the absorbers (see previous 
section): a\ x cm sa a 2 neutral + a 2 clouds . 

Of course, a significant fraction of the gas whose turbu- 
lent velocity is measured by radio observations is behind 
the ionizing stars in front of which absorption lines from 
neutral elements originate, hence the velocity structure in 
our case could be different (Pettini & Lipman 1995). The 
effect of the determination of the turbulent velocity on 
column densities is discussed in Sect. I4.1.TH 
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Table 3. Column densities derived in IZw 36. Errors are given 
at 2a, upper limits at 3a. Upper limits for H2 are calculated for 
the same turbulent velocity as the neutral phase (12.2 km s _1 ). 
The Fe ill column density is calculated from the Fe ill A1122.53 
line which could originate from both the ISM and stellar at- 
mospheres. 



Species 


log iV(cni _2 ) 


log (X/H) 


[X/H] a 


Hi 


21.30+^ 






Cii 


is.3i±g : iS 






Cm 


14.38i«;i 8 9 






Ni 


14.42+°;! 


-6.88«; 3 3 ? 


71 +0.47 
z - ' 1 -0.42 


Nil 


14.02lJ;SS 






Oi 


16.80+1 g 


-4.50±l;?g 


-1.191?;?° 


Ovi 


13.97i«; 2 6 2 






Sin 


15.34t°;H 


-5-96±g : g? 


i - du -0.39 


Pn 


13.00+. ,; 2 * 


-8.3o±°:SS 


-1 76+ ' 41 

L - '"-0.51 


Sm 


14 84+0- 4 3 






An 


< 13.30 


<-8.00 


<-2.55 


Fen 


14.511°;^ 


-6 79+ 031 

°- ' y -0.23 


-2 26 +0 - 34 

^• zu -0.26 


Fein 


14.47±g;il 







H 
H 
H 



2 (J=0) 
2 (J=l) 
1J=2)_ 



< 14.41 

< 14.61 

< 14.32 



[X/Y]=log (X/Y)- 
solar value. We 1 
Lodders (2003). 



log (X/Y) , where log (X/Y) Q is the 
se the recommended solar values of 



Table 4. Effects of the total line broadening atot (in pixels, 
see Sect. l3.lt on column densities (log iV(cm~ 2 )). Note that if 
a tot is overestimated so are the column densities. 







<T iot = 12 


15 


21 


24 


27 


log 


iV(Ni) 


14.30 


14.30 


14.32 


14.36 


14.42 


log 


AT(Siii) 


15.08 


15.08 


15.15 


15.20 


15.34 


log 


JV(Pli) 


12.99 


13.02 


13.14 


13.16 


13.17 


log 


AT (Feu) 


14.20 


14.20 


14.26 


14.30 


14.49 



tested several values of atot in order to test its impact on 
the Fe II, Si II, P II, and N 1 column densities. It can be 
seen in Table 0] that for atot between 12 and 27 pixels, 
column densities do not vary by more than a factor 2.5 
in the worst case (Si 11) . A bad determination of the line 
broadening would not change the results of Sect. 03 

4.1.1. Molecular hydrogen 

No H 2 lines were detected at the redshift of IZw 36 
(Table EJ, as for IZw 18 (Vidal-Madjar et al. 2000), 
Markarian 59 (Thuan et al. 2002), and SBS 0335-052 
(Thuan et al. 2003). There are several possible physical 
explanations for this absence as pointed out by Vidal- 
Madjar et al. (2000): the ionizing flux of the massive 
stars which can be high enough to destroy H2 molecules, 




8x10 - 



Fig. 1. Hi and Fen lines in IZw 36. The flux is in 
erg s _1 cm -2 A -1 and the observed wavelength in A. The data 
(histogram) are binned by a factor 3 for display purposes. The 
thin line represents the line profile before the convolution with 
the total line broadening represented by atot in Sect. I3.l1 and 
the thick line is the convolved profile. For each of the spectral 
windows, we used the two observations and all the available 
FUSE detection channels (Sect.|5J to perform the fit. Here we 
show the spectrum of segment LiFl for the upper plot and 
of segment LiF2 for the lower plot, a) upper Plot of the Ly/3 
line. The center of the line is blended with terrestrial airglow 
emission lines. Atomic hydrogen in the Milky Way is invisi- 
ble because of its low column density ~10 19 cm -2 (estimated 
from Heiles 1975). The Sin A1020.70 lines and Ovi A1031.93 
lines from the Milky Way and IZw 36 are also visible in this 
spectrum, b) lower Plot of the strong Fell Al 144.94 line. 



the scarcity of dust grains on which H2 molecules are 
synthetized, and a low Hi volumic density. Also, Vidal- 
Madjar et al. (2000) argue that the molecular medium 
could be clumpy so that the few clumps in the lines of sight 
of the many blue stars in the entrance aperture would not 
be detected. 

4.1.2. Hi 

The Hi column density has been calculated using all the 
available Lyman lines (Ly/3, Lye, Ly£, Ly?7, and Ly(9) but 
is mainly constrained by the damping wings of the Ly/3 
line (Fig. ^,) since the other H 1 lines are strongly satu- 
rated. We used a 30 A wide window centered on Ly/3 to 
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constrain the continuum by a 4-order polynomial. Wc find 
log iV(Hl)= 2l.30t° Q ° 1 9 Q - The errors on log AT(Hi) are the 
smallest amongst the specie we investigate. This is due 
to the damping profile of the Ly/3 line which is observed 
through both LiF and SiC channels. 

Our value is somewhat higher than log iV(Hi)=21.15 
found by V83 who used only the red damping wing of the 
Lya line from IUE observations, but is certainly better 
constrained in our case because of the two damping wings 
of the Ly/J line and because of the better resolution and 
the better signal-to-noise ratio of our observations. 

4.1.3. Ni and Nil 

N I is the most useful species in our study since most 
of it must be present in the neutral ISM. Six Ni lines 
were used (Table These lines are not saturated so 
that the determination of log 7V(N i) is not strongly de- 
pendent on the turbulent velocity b (Fig. [2J. We obtain 
log iV(Nl)=14.42l°:|t 




5 10 15 20 25 

b(lms-l) 



Fig. 2. Effects of the determination of the turbulent veloc- 
ity b on the Ni column density (solid line). Errors in b do 
not affect much the column density determination. The A\ 2 
curve (dashed line) is relatively narrow and is minimal for 
6=12.2 km s~\ At 2a (A X 2 =4) we obtain b=12.2±|* corre- 
sponding to log N(N i)=14.43±g;o|. 

The A1083.99 Nil line is barely detected since it 
is located in a spectral region which is very noisy. 
We are however able to calculate a column den- 
sity: log jV(Nn) = 14.02±J;£|. 

4.1.4. Oi and Pn 

The three O I lines in IZw 36 are either blended or strongly 
saturated (Table |2J, resulting in large errors (Table EJ). 

The A963.83 Pn line has a large oscillator strength 
(/=1.25) and is detectable even for low Pn column den- 
sities. In IZw 36 spectra, this line is blended with N I lines 
from the Milky Way and from IZw 36. Given the fact that 
the Ni column density is well determined and that the 
Pn line is not saturated, the blending can be corrected 



All 




1047 1048 1049 1050 1051 

wavelength 



Fig. 3. Lack of argon in IZw 36. See Fig. Q for a descrip- 
tion of the plot. Data are binned by a factor 2 for display 
purposes. The line plotted at the redshift of IZw 36 is for 
log AT(Arl)=14.0 and 6=12.0 km s _1 . 



for. We also use the A1152.82 Pn line with a fainter os- 
cillator strength to constrain the final P n column density 
as log iV(Pn)= 13.00^;^. 

4.1.5. Sin, An, Fen 

The Sin column density is only given by the A 1020. 70 
line which is slighly saturated. Unfortunately, the other 
strong Sill line at 989.90 A is contaminated by terrestrial 
airglows. 

The other species lead to reliable column densities: Ar I 
(Fig. El for which we have a good upper limit, and Fen 
(Fig. ^>) with a large collection of lines. 

4.2. Abundances in the neutral region 

In order to derive total abundances from ionic column 
densities, it is necessary to know the ionization states of 
each species in the neutral region. 

In principle, in the Hi region we expect to find all ele- 
ments with an ionization potential larger than that of hy- 
drogen (13.6 eV) as neutral atoms. Elements with smaller 
ionization potentials must be found as single-charged ions. 
For example, we can safely assume that iron and silicon 
exist mainly as Fe II and Si n in the neutral region and that 
the Fen and Sin lines are essentially produced there. We 
have checked that the column densities of these ions are 
considerably smaller in the Hil region using the model 
E2E1 from the grid of Stasihska (1990). The model E2E1 
agrees the best with the physical parameters of the Hil 
region derived by V83. The resulting ionic fractions of Si n 
and Fen in the Hil region in the best fitting model are 
only a few percents. 

The situation is however not so simple for ele- 
ments with an ionization potential not much larger than 
that of hydrogen, because UV photons able to pro- 
duce their ionization may be present in low-density, 
partly-ionized regions of the ISM. For example, An 
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(ionization potential 15.8 eV) is easy to ionize be- 
cause of its relatively large photoionization cross section 
(Sofia & Jenkins 1998). Therefore, in regions with an ion- 
ization degree n e /n#«0.2, the actual Ar/H can be larger 
by 0.2 to 0.7 dex than An/Hi, depending on the hard- 
ness of the far-UV radiation. The situation is similar but 
less severe for Oi and Ni (ionization potentials 13.6 eV 
and 14.5 eV respectively). Indeed, in the conditions of 
the Local Cloud for which the correction for Arl/Hl is 
0.36 dex, the corresponding correction for Ni/Hi is only 
0.05 dex (Sofia & Jenkins 1998). However the ionization 
correction in BCDs could be different and could reach 0.2- 
0.3 dex. We will first neglect this correction, and because 
little N I and O I are expected in the H n region, we assume 
that the Ni and Oi lines are representative of the total 
nitrogen and oxygen abundances in the neutral region. 

The abundances of the neutral region can now be com- 
pared with the ionized region values. For the comparison, 
we use metal abundances relative to hydrogen (X/H). 

5. Comparison with the ionized gas 

Figure 2] compares several abundances ratios representa- 
tive of the ionized phase of IZw 36, derived from optical 
and UV emission lines, with those of the neutral gas. 
We define the underabundance Sni(x), as the logarithmic 
difference of an abundance ratio x between the H n and 
the H I regions by : 

S m (x) = log {x)mi - log (x)m 



If 

-5- -r - 

s ,-- I I - 

- 8 i i 

: N/H O/H Si/H Ar/H Fe/H : 
-9 C J 

Fig. 4. Comparison of X/H between Hi and Hn regions 
in IZw 36. Error bars are at 2a. Squares: this study (as- 
suming the abundances of Sect. 14.21 and using An for ar- 
gon i.e. without any ionization correction). Triangles: V83. 
Circles: Izotov et al. (1997) for N, O, Ar, and Fe. The sili- 
con abundance in the ionized region, not measured in IZw 36 
is assumed to be equal to — 1.5±0.1 which is the mean value 
observed in BCDs (Izotov & Thuan 1999). 



Si II column density is uncertain. Within the large error 
bars, the abundances in the ionized region and in the neu- 
tral region are in agreement. 

Because of all these uncertainties, and because silicon 
and iron can be depleted on dust grains, we prefer to focus 
on nitrogen abundance which is much more reliable. 



The column densities of oxygen and silicon of IZw 36 
are similar in the neutral and in the ionized gas within 
large uncertainties. However the result for nitrogen is 
much more reliable, and is certainly the best established 
result of this study, yielding <5hi(N/H)=1.2±0.3; in other 
words N is about fifteen times less abundant in the neu- 
tral gas with respect to the ionized gas. This differ- 
ence is larger than 3o\ Other metals for which abun- 
dances are rather accurately determined also show under- 
abundances: <5 H i(Fe/H)=0.7±0.2 and <5 m (Ar/H)>1.5±0.1 
(without any attempt to correct for partly-ionized, low- 
density regions). 

The following discussion aims to identify the effects re- 
sponsible for the observed iron, nitrogen, argon, and oxy- 
gen underabundances. 

5.1. Iron and silicon 

A concern with the abundance of iron in Hll regions is 
that it is derived from the weak A4658 line of Fein, and 
is relatively uncertain. In any case, we find that iron is 
underabundant by <5 H i(Fe/H)=0.7±0.2 in IZw 36. This 
underabundance might be overestimated since it is likely 
that iron is more efficiently depleted on dust grains in the 
neutral phase than in the ionized phase. 

In IZw 36, the only one available Sin line is slightly 
saturated (Sect. I4.1.5fl so that the determination of the 



5.2. Nitrogen 

The underabundance in IZw 36 #hi(N/H)=1.2±0.3 is too 
large to be accounted for by chemical inhomogeneities 
alone in the neutral gas, given the underabundance of iron. 
In order to illustrate the underabundance of nitrogen in 
the neutral gas, we have compared in Fig. our best fit 
of the A1134.98 Ni line, one of the 6 lines we used, with 
the profile assuming the N/H ratio as in the ionized gas. 
The two profiles are inconsistent within 3<r. 

Furthermore, even including the largest possible ion- 
ization correction (i.e. including N n and taking the upper 
error bars on N I and N II column densities) , the under- 
abundance of nitrogen is still positive: fer(N/H)=1.0±^. 
This supposes that 75% of the nitrogen is ionized into N n 
in the neutral gas of IZw 36. Although the determination 
of Nil column density is uncertain (Sect. I4.1.3fl . we can 
reasonnably conclude that nitrogen is indeed underabun- 
dant in the H I gas of IZw 36. 

5.3. Argon 

The argon deficiency is ($Hi(Ar/H)>1.5±0.1. Argon is pro- 
duced in massive stars and should not be depleted in 
a low-density Hi cloud (Jenkins et al. 2000). A major 
part of the underabundance in the neutral gas is certainly 
due to an ionization correction from low-density, partly- 
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Fig. 5. Fit of the A1134.98 line. Data are binned by 
a factor 2 for display purposes. The solid line is for 
log (N i/H i)=— 6.88 and shows our best fit. The dashed line 
is for log (N i/H i)=— 5.63 and does not fit well the data. 



ionized regions of the ISM (Sect. I4.2|l such as argon is 
mainly ionized into Aril. 

5.4. Oxygen 

The O /H ratio in the neutral region of IZw 36 is difficult 
to assess because of the large errors on the O I column 
density (Table EJ). However, it is possible to estimate its 
value using other elements which are produced in the same 
massive stars. Phosphorus turns out to be a good tracer of 
oxygen (Lebouteiller 2003), since argon is partly-ionized 
in the neutral region (Sect. 15.3(1 and since the depletion 
of silicon on dust grains is uncertain. Depletion is not a 
concern for phosphorus since it is not much depleted, and 
in the same way as oxygen (Savage & Sembach 1996). 

We calculate O p which estimates the neutral 
oxygen abundance assuming [P/O]~0 in the neu- 
tral gas (using the recommended solar value of 
Loddcrs (2003) log (P/O) =-3.23±O.O9). We find 
log (O p /H)«— 5.1±0.5, giving an underabundance 
<5m(O p /H)«0.9±0.5. Within the errors, this estimation 
of the oxygen abundance is lower with respect to that in 
the ionized gas. Given the low depletion of phosphorus 
and oxygen in the Galactic ISM (Andre et al. 2003 and 
Savage & Sembach 1996), it is very likely that oxygen is 
genuinely deficient in the neutral medium of IZw 36. 

Assuming log (O p /H)=— 5.1±0.5 in the Hi gas, 
we find log (N/O p )=— 1.8±0.6 (2a uncertainty) lower 
than, although consistent with, the value — 1.49±0.01 
in the ionized gas (Izotov & Thuan 1999). Also, we 
find log (Ar/O p )<— 3.0 as compared to the solar value 
log (Ar/O) =-2.1±O.l (Lodders 2003), confirming that 
argon is likely ionized into Ar II as stressed in Sect. 15.31 



6. Conclusions 

Several heavy elements in the neutral gas of IZw 36 are un- 
derabundant in the neutral gas with respect to the ionized 
gas. Various effects can be responsible for the underabun- 
dances <5hi(X/H): 

— The presence of an unprocessed neutral gas, less chem- 
ically evolved, in the line of sight can reduce all the abun- 
dances X/H in the neutral gas with respect to the ionized 
gas. 

— A more efficient depletion on dust grains in the neutral 
medium can account for the iron deficiency. 

— An additional metal enrichment could be responsible 
for nitrogen and oxygen (using phosphorus as a tracer) 
overabundance in the ionized gas, although in this case, 
silicon which is an a-element should be also overabundant 
(but see caption of Fig [J}. 

Our findings that metals are observed in the neutral 
region of IZw 36 indicates that this gas phase has been al- 
ready enriched by previous star formation, either quiescent 
or in bursts (Legrand 2000). Schulte-Ladbeck et al. (2001) 
detect indeed an old stellar population which could have 
enriched the neutral region of IZw 36. 

In this study, we considered Hn region abundances 
derived from forbidden lines in the visible and the ultravi- 
olet. However, the abundances derived in this way might 
be underestimated by factors (>2 for oxygen for instance) 
as claimed by Tsamis et al. (2003) compared to the more 
reliable abundances deduced from optical recombination 
lines or far-IR forbidden lines (which unfortunately are not 
available for the BCDs). Consequently, the abundances 
differences between the two gas phases of IZw 36 could be 
even larger than those given in Fig. 0] 

At present time, the overall picture still remains un- 
clear. To assess which of the effects (unprocessed neutral 
gas, depletion on dust grains, metal enrichment) are dom- 
inant, it is necessary to investigate the neutral gas of more 
gas-rich star-forming galaxies. Furthermore, with a larger 
sample, we can reasonnably expect to find several good 
oxygen determinations. 
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